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The most common intrinsic defects in dislocation-free float zone (FZ) silicon 
crystals are the A- and B-type "swirl** defects. The mechanisms of their 
formation and annihilation have been extensively studied by Petroff and 
de Rock^D , Foell, Goesele and Kolbesen^) t and others. Another type of 
defect in dislocation-free FZ crystals, found by Roksnoer and van den Boom*-*), 
is referred to as a D-type defect. Concentrations of these defects can be 
minimized by optimizing the growth conditions and the residual swirls can be 
reduced by the post-growth extrinsic gettering process. 

Czochralski (Cz) silicon wafers are known to exhibit higher resistance to slip 
and warpage due to thermal stress than do FZ wafers. Sumino and Yonenaga^ 4 ' 
examined various conditions under which dislocation-free Cz crystals have 
higher mechanical strength than dislocation-free FZ crystals. They observed 
no difference in the yield strength when external damage had been removed by 
chemical etching. The Cz crystals containing dislocations are more resistant 
to dislocation movement than dislocated FZ crystals because of the locking of 
dislocations by oxygen atoms present in the Cz crystals. Nitrogen doping of 
FZ wafers at concentrations as low as 1.5 x lO* 5 cm - -* was found by Sumino, et 
al.* 5 ) to make such wafers less susceptible to slip and warpage than undoped 
Cz and FZ wafers containing dislocations. 

Recently we have applied a transverse magnetic field during the FZ growth of 
extrinsic silicon.* The objective was the study of the spatial distribution 
of the major dopant, under the assumption that the magnetic field would reduce 
fluctuations in the flow occurring in the melt and thus in the dopant 
concentration. Important changes in spatial distribution of the dopant were 
indeed observed, but the level of fluctuations was not substantially 
decreased. The observed flow patterns, as revealed by striation etching and 
spreading resistance in Ga-doped silicon crystals, indicate strong effects of 
the transverse magnetic field on the circulation within the melt. At fields 
of 5500 gauss, the fluid flow in the melt volume is so altered as to affect 
the morphology of the growing crystal. For crystals grown without rotation, 
the cross section of the crystal becomes elliptical, with the major axis of 
the ellipse aligned along the field direction. Melt flow in the directions 
parallel and perpendicular to the field are distinctly different. We believe 
that the magnetic field offers a unique experimental tool to help elucidate 
the mass and energy transport occurring in the float zone process. 

* This work has been supported in part by the U.S. Air Force Materials 
Laboratory. 
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Introduction 


• FLOAT ZONE (FZ) Si CRYSTAL GROWTH 

• CRYSTALLINE DEFECTS 

• DEFECT GETTERING 

• MECHANICAL STRENGTH OF FZ CRYSTALS 

Growth of <100> Dislocation-Free FZ Si Crystals 

PARAMETERS AFFECTING THE DISTF;!BUTION AND CONCENTRATION 
OF INTRINSIC DEFECTS 

• PULL RATE 

• ROTATION RATE OF CRYSTAL 

• AMBIENT ATMOSPHERE 

• COOLING RATE OF CRYSTAL 

• MAGNETIC FIELD 

Observable Imperfections in Dislocation-Free 
FZ Si Crystals 

• STRI ATIONS 

• "SWIRLS" 

• STACKING FAULTS 
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Swirls 


DENSELY PACKED DISCRETE DEFECTS 

A-, B-, AND D-TYPE DEFECTS IN FZ DISLOCATION-FREE CRYSTALS 

A-TYPE 

• DISLOCATION LOOPS OR CLUSTERS OF DISLOCATION LOOPS WHICH 
ARE FORMED BY AGGLOMERATION OF Si SELF-INTERSTITIALS 

• DISTRIBUTED IN STRIATED PATTERN NEAR BANDS OF IMPURITIES 
WITH DISTRIBUTION COEFFICIENT LESS THAN ONE 

• SIZE: ~1 - 40^m 

• CONCENTRATION: 10 7 - 10 9 cm -3 (H. FOLL, et al.) 

• PULL RATES OF CRYSTAL FOR FORMATION OF DEFECTS: 0.2 - 
4.5 mm/min (23 mm diam) 

B-TYPE 

• THREE DIMENSIONAL AGGLOMERATES OF Si SELF-INTERSTITIALS 
AND IMPURITY ATOMS, PARTICULARLY C 

• DISTRIBUTED IN STRIATED PATTERN SIMILARLY TO A-TYPE SWIRLS 

• SIZE: < 30 nm 

• CONCENTRATION: <~ 10 11 cm -3 (H. FOLL, et al.) 

• PULL RATES OF CRYSTAL FOR FORMATION OF DEFECTS: 0.2 - 
5 mm/min (23 mm diam) 

D-TYPE 

• CONDENSATION OF POINT DEFECTS (VACANCIES), WHICH DO 
NOT PRECIPITATE ON CARBON NUCLEI 

• HOMOGENEOUSLY DISTRIBUTED IN THE CRYSTAL 

• SIZE: 30 nm - 1 /xm 

• CONCENTRATION: 10 9 cm -3 (P. J. ROKSNOER, et al.) 

• PULL RATES OF CRYSTAL FOR FORMATION OF DEFECTS: 

>6 mm/min 
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CONCENTRATION, cm 


Swirl Concentration vs Growth Rate 



FROM: P. J. ROKSNOER AND M. M. B. van den BOOM 

J. OF CRYSTAL GROWTH 53 (1981) 563 - 573 


Extrinsic Gettering of Swirl Defects 

ABRASION GETTERING (C. L. REED AND K. M, MAR) 

• MATERIAL 

<11 1>, P-DOPED, FZ Si CRYSTAL WITH SWIRL DEFECTS 

WAFERS 305±13 jxm THICK 

SPECULAR DAMAGE-FREE FRONT SURFACE 
CHEMICALLY ETCHED STRAIN-FREE BACK SURFACE 

• ABRASION 

BACK SIDE DAMAGE APPLIED BY MECHANICALLY LAPPING WITH 
2-15 fim SIZE DIAMOND AND ALUMINA PARTICLES 

• OXIDATION AND CHEMICAL ETCHING 

OXIDIZED AT 1150°C TO ACTIVATE DEFECTS AND WRIGHT ETCHED 
TO REVEAL DEFECTS 

• RESULTS 

ABRASION GETTERING REDUCED THE ETCH PIT DENSITY FROM 
5 x 10 5 TO< 5 x 10 2 cm' 2 
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Mechanical Strength of FZ and Cz Si Crystals 
(K. Sumino et al, 1981) 


• YIELD STRENGTH SAME FOR DISLOCATION-FREE CHEMICALLY 
POLISHED FZ AND CZ WAFERS 


• DIFFERENCE IN MECHANICAL STRENGTH DUE TO DISLOCATION 
PRIOR TO DEFORMATION 


• DISLOCATIONS NOT GENERATED IN CZ FOR STRESSES T<T C 


• DISLOCATIONS GENERATED IN FZ FOR VERY LOW T 


Effect of Nitrogen on Mechanical Strength 
of FZ Crystals 
(K. Sumino et al, 1 983) 

• 1.5 x 10 15 cm" 3 N-DOPED FZ LESS SUSCEPTIBLE TO SLIP 


• N ELECTRICALLY INACTIVE IN Si. OCCUPIES INTERSTITIAL SITE 

• DISLOCATION-FREE FZ, CZ HAVE SIMILAR STRESS-STRAIN 
CHARACTERISTICS 

• DISLOCATED FZ N-DOPED HAVE HIGH UPPER YIELD STRESS 

• DISLOCATIONS IN FZ N-DOPED IMMOBILIZED UNDER LOW STRESSES 
AT ELEVATED TEMPERATURES 
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Introduction 


• EXPERIMENTAL APPARATUS FOR APPLYING TRANSVERSE 
MAGNETIC FIELD TO FZ MOLTEN ZONE 

• TYPICAL RESULTS FROM ROTATING CRYSTALS 

• UNUSUAL RESULTS FOR NON-ROTATING CRYSTALS 

• SIGNIFANCE OF RESULTS 


Experimental Parameters 


• CRYSTAL ORIENTATION 

• CRYSTAL DIAMETER 

• PULL RATES 

• ROTATIONAL RATES 

• MAGNETIC FIELD STRENGTH 


<100> AND <11 1 > 
UP TO 30 mm 
2 TO 4 mm/min 
OTO 12 rpm 
0 TO 5500 G 
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Striation Etch Pattern, Crystal Z309; 5000 Gauss 
Transverse Magnetic Field (12 rpm Section) 





n, ROOT cm/RAVG 


Typical Spreading Resistance Data 



0.10 0.20 0.30 0.40 0.50 0.60 


Power Spectral Densities 
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ORIGINAL PAGE It 
OF POOR QUALITY 

5500-Gauss Field Applied to Rotating Crystal 

OBSERVATION HYPOTHESIZED CAUSE 


• GROWTH INTERFACE "FLATTENS" 

• CRYSTAL DIAMETER INCREASES 

• DOPANT CONCENTRATION INCREASES 

• ROTATIONAL STRIAE ENHANCED 

• FINE STRIAE ALTERED 


ALTERED HEAT FLOW PATTERNS 

FLATTENING OF INTERFACE 

THICKENED DIFFUSION 
BOUNDARY LAYER 

REDUCED FLUCTUATIONS 
REDUCED "TURBULENCE" 


Crystal Grown With No Rotation in 
5000 Gauss Transverse Field (Z321) 



FIELD DIRECTION 

A 




A 

FIELD DIRECTION 
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5 500-Gauss Field Applied to Non-Rotating Crystal 


OBSERVATION 


HYPOTHESIZED CAUSE 


• DOPANT CONC INCREASES 

• ROTATIONAL STRIAE ABSENT 

• FINE STRIAE ALTERED 

• SR FLUCTUATIONS REDUCED 
FOR F <300 cm" 1 (~2 Hz) 

• CROSS SECTION — ► ELLIPTICAL 

• GROWTH INTERFACE— •— CYLINDRICAL 


THICKENED DIFFUSION 
BOUNDARY LAYER 

AS EXPECTED FOR 0 RPM 

REASON NOT CLEAR YET 

INCREASED MELT "VISCOSITY" 

Bx V EFFECTS 
ALTERED HEAT TRANSFER 


Summary 

• MAGNETIC FIELD IMPACTS CIRCULATION PATTERNS IN THE 
MOLTEN ZONE 


• FORCES FROM FIELD ARE LARGE ENOUGH TO PRODUCE 
USEFUL EFFECTS 


• FOURIER ANALYSIS OF SPREADING RESISTANCE MAY REVEAL 
FLUCTUATION SPECTRA IN MELT 


• MAGNETIC FIELD PROVIDES NEW TOOL TO STUDY THE FZ 
GROWTH PROCESS 
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DISCUSSION 


NEITZEL: 

ROBERTSON: 

NEITZEL: 

ROBERTSON : 


NEITZEL: 


I've heard the word turbulence mentioned a couple of times today and 
I think we have to be very careful when we mention the word 
turbulence because it implies something definite, in a fluid dynamic 
sense, that I don’t think is occurring in your geometry for your 
region of parameter space. 

Do you have some good reason for saying it's not occurring? 

I would guess it's not. We did a rough calculation for the maximum 
rotational Reynolds number that you should see in your 
configuration. It's probably on the order of 1000. For the 
transition to turbulence, the transitional Reynolds number is 
3 x 10 5 . 

All I can tell you is that the flows in that melt volume are very 
strong. They are being driven by electromagnetic forces. They are 
not being driven by convective forces or Marangoni surface-tension 
forces or the rotational forces. Those are all there, but 
superimposed on that is a circulation pattern from the 
electrodynamic forces. Professor Muhlbauer has just published 
information that indicates that that force is 50 to 100 times 
stronger than these other forces. If you look at particles in that 
melt zone, for example, they are in violent motion. 

I don't doubt that, but I just want to say that it's possible to have 
very oscillatory laminar flows that are not turbulent flows. 
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